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BODIES OF RE70LUTI0R AT HIGH SUPERSOHIC AIRSPEEDS 


By A. J. Eggers, Jr., and Raymond C. Savin 


SUMMARY 


Flow at Mgji supersonic speeds about a body of revolution is 
investigated analytically. With the assumption that the flow at the 
vertex is conical, it is found that algebraic solutions can be obtained 
which yield the Mach numbers and pressures at the surface over a con- 
siderable range of free-stream Mach numbers and apex angles. In the 
special case of cones, these solutions define the entire flow field with 
good accuracy, aufl may therefore provide a useful adjunct to the well- 
knovm M.I.T. tables. 

The investigation of flow downstream of the vertex reveals that 
when the value of the hypersonic similarity parameter for the flow 
(i.e., the ratio of the free-stream Mach number to the slenderness ratio 
of the body) is large con 5 )ared to 1, the Mach number along a streamline 
(downstream, of the nose shock) varies with flow inclination angle in 
approximately the same manner as for two-dimensional (Prandtl -Meyer) 
flow. In the special case of streamlines near the surface, it is si 2 g- 
gested that this parameter may approach 1. This resTolt an d the solutions 
obtained for flow at -the vertex are combined to yield ■»diat might be 
cal led a conical-shock -expansion method for calculating the Mach number 
and pressure distributions at the surface of a body. These calculations 
are shown to be particularly simple in the case of slender bodies. 

Surface Meich number and pressure distributions calculated with the 
simplified methods of this paper for a number of ogives are fo-und to be 
in good agreement with those obtained With the method of characteristics 
. at valiies of the hypersonic similarity parameter greater than 1. In the 
case of the conical-shock -expans ion calculations, the agreement is within 
the oarder of accuracy of the characteristics solutions When the hyper- 
sonic similarity parameter has a value of only 2. Because of the rela- 
tive simplicity of these ceilculations, the methods for determining the 
flow at the siirface of a body mEiy prove useful for engineering purposes. 
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INTRODUCTION 


Determination of supersonic flow fields about nonlifting bodies of 
revolution by mea n s of tbe method of characteristics (see, e.g., refer- 
ence l) is generally accepted as an accurate but tedious and time- 
consimiing operation. Because of the latter features of this method, 
recourse is often made to sinpler methods ■which, altho\agh less accurate, 
can be applied with relative simplicity and rapidity. 

Perhaps the most ■widely used sinplified theory of axieilly symmetric 
s^upersonic flow is the linear theory, ■sdiich ■vra.s first employed by Kdrmdn 
and Moore (reference 2) to study sipersonic flow about a body of revolu- 
tion. This theory has a limited rainge of applicability, ho^wever, due to 
the assunptlon in its development of potential flow ■with infinitesimal 
distinrbances . In particular, it is only applicable ■with good accuracy 
to bodies of practical slenderness ratios operating at low supersonic 
airspeeds . 

More recently. Van Dyke (reference 3) has developed a second-order 
s^upersonic flow theory which srields results "which are generally more 
acc^u 2 *ate than those ob"tained ■with the linear or first-order theory, 
altho\:igh the calcxilations are, of course, more lengthy. The range of 
applicabili"ty of Van Dyke’s second -order solution is specifically 
limited, however, to cases for ■which the product of the free -stream 
Mach number and maxinnmi slope of the body is less than 1; hence, it is 
not generally suitable for calcialating the flow about bodies of usual 
proportions at high stpersonic, or hypersonic, airspeeds . 

Flow at hypersonic airspeeds in the limit as the free -stream Mach 
n\miber approaches infinity eind the ratio of specific heats of the gas 
do"wnstream of a shock "wave approaches 1 has been st'udied by Busemann 
(reference 4) . It was found that under these circums^tances a very simple 
expression is obtained for the pressures acting on a nonlifting body of 
revolution. This e3q)ression provides in general, however, only qualita- 
tively acc^urate press\ire distrlbxitions at high but finite free -stream 
Jiach numbers ■whei^ the ratio of specific heats is greater than 1 and is 
usu all y closer to the ideal diatomic gas value of 1.4. 

Perhaps the first step in the direction toward pro"vldlng simplified 
methods for determining axially symmetric flows in the high s^upersonic 
speed range ■was made by Tsien (reference 5) • Tsien demonstrated that a 
similarity law (■which is analogo'us "to the well-known law applicable at 
low s^upersonic speeds) exists for slender pointed bodies in this range. 
Thus, "With the aid of this law, the flow about a family of shapes a-1 1 
having the same thickness distribution can be easily determined, pro"vlded 
the ratio of the free -stream Mach niimber to the slenderness ratio is the 
same for all shapes, and provided the flow about one of these shapes is 
known. This result materially reduces the net labor associa"ted ■with 
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calculating flows atout a number of related bodies of revolution opera- 
ting at high supersonic airspeedsj bowever, tbe problem remains of 
calculating tbe flow about a representative body for each value of tbe 
similarity parameter. At present, tbis problem can only be solved with 
good accuracy by means of tbe method of ebaracteristics . 

There appears, then, to be a need for a sinplified theory which 
provides, with engineering accuracy, tbe Ifecb number and pressure distri- 
butions about pointed nonlifting bodies of revolution operating at high 
supersonic airspeeds. Tbe development of such a theory is undertaken 
in tbis report, and it is found that by treating tbe flow in two parts, 
namely, tbe flow at tbe vertex and tbe flow downstream of tbe vertex, a 
theory of tbe desired sinplicity and accuracy can, in fact, be obtained. 
In tbe special case of flow about cones, tbe theory demonstrates sur- 
prising accuis,cy over a considerable range of stpersondc airspeeds. 
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SYMBOLS 


first family characteristic line 
second family characteristic line 
maximum diameter of a body of revolution 
entropy 

hypersonic similarity parameter 

characteristic body length (measured from vertex to most foiward 
point of maximum diameter) 

Mach nuinber (ratio of locetl velocity to local velocity of sound) 

distance measured normal to ray passing tbrou^ vertex of cone 

static pressure 

dynam i c pressure ' 

pressure coefficient { ? . ~?0 

\ ^ 




gas constant 

radial distance from vertex of tbe cone 
distance along a streamline 
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u,v 

V 

/S. 

V 

P 

7 

6 angle of flow Inclination with respect to the body axis 
p mass density 

01 angle between axis of cone and. ray passing through vertex of cone 

Subscripts 

o free -stream conditions 

c cond.it ions on the sttrface of a cone 

IT conditions on the siirface at the vertex of a body 

s cond.itions Immediately behind the shock wave at the vertex of a 

body 

AMLYBIS OF FLOW ABOUT A NONLIFTHTG BODY OF EEVOLUTIOW 


velocity con^ionents parallel and normal, respectively, to ray 
passing through vertex of cone 


resultant velocity ^/u^ + 

TnaYlrniim velocity obtainable by expanding to zero tenperature 
rectangular coordinates 


.(arc sin 


Mach angle i arc sin — 

' M 

y 

ratio of specific heat at constant pressure to specific heat at 
constant volume 


This investigation is concenied with a simplified method of calcu- 
lating axially symmetric flow about a body of revolution traveling at 
high supersonic airspeeds. It is assumed throughout the analysis that 
the disturbed flow is everywhere supersonic, and thus, of course, that 
the body has a sharp nose or vertex. With these restrictions on the 
free -stream Mach number and body shape, it is evident that a typical 
flow. field will be characterized by the bow shock wave lying close to 
the surface of the body. For the purpose of analysis, it is convenient 
to study the flow field in two parts, namely, the flow at the vertex 
and the flow downstream of the vertex. The results of these two phases 
of the investigation will then be applied .to the determination of the 
Mach number and pressure distributions over the body surface. 
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Flow at the Vertex of a Body 


It follows from the assumptions basic to this analysis that the 
nose shock is attached., and. consequently that the initial conditions of 
the flow at the vertex will he the same as for a cone tangent to the 
body at the vertex and. operating in the same free stream. In the classi- 
cal paper of Taylor and. MaccoU (reference 6) the basic differential 
eqmtions were developed, for steady, axially symmetric supersonic flow 
about cones. ^ The essential assunption in their analysis was that fluid, 
properties were constant along radial lines passing through the vertex; 
hence the flow was irrotational. This assuii 5 )tion will be en 5 >loyed. here 
to red.evelop these differential equations in a form more suitable for 
this analysis. 

A schematic diagram of axially symmetric supersonic flow about a 
cone is shown in figuire 1. Considering first the flow downstream of the 
shock wave, the equations of motion and continuity are written (after 
reference 6) in the fonns 

( 1 ) 

and 

— (p V sin to) + 2 p u sin w=0 (2) 


(3) 


and 

V = -V sin (u-8) (5) 

Therefore, equations (l) and (2) may be combined to obtain the relation 


As pointed out in reference 6, Bxisemann had previously suggested a 
graphical solution of the problem of supersonic flow past a cone. 


respectively. The condition of irrotationality is given by the 
expression 



Now it is evident from figure 1 that 

u = V cos (to-6) 
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2^ (v2 - y2) V cos (w-8) + ^^-sin (o)-8) 


Y- coto) + 


do ) J 


+V^-^ sin (o)-8) = 0 
dw ^ 


( 6 ) 


Similarly, tlie irrotationality et^uation may le "written 


dV 


- "tan (ci)-6) dS 


(7) 


Comibining equations (6) and (7), there results, "Haen, the equation of 
motion 


Zri 

2 


" 'VV T ^ 
Y I 


cot (‘*1-8) ^1 + -cot (0 + tan 


(o>-6) ^ 
do) 


- "ban (o)-8) ~ = 0 
dco 


(8) 


"Which "will he employed in the following development. 

Equations (7) and (8) are amenable to nimerical integration follow"- 
ing methods similar to those emplos^ed earlier in reference 6. In order 
to obtain algebraic solutions to flown at the ver"tex, however, it is 
necessary to stnqjlify these eqmtions. To this end, since high Mach 
nimiber flows are of principal in"fcerest, an assun^ition "which has proved 
useful in s"tudying other aspects of these flows (see, e.g., reference 7) 
is employed, namely, the magni"tude of the reB"ul"bant veloci"fcy vector is 
essentially cons"fcant, the chaJige in veloci"fcy being primarily one of 
change in direction. With the aid of this asstimption it is possible to 
first detezTnine a relation be"tween 5 and u for a given V using 
equation (8), and then to determine more accura"tely "the nature of the 
dependence of V on 6 and u using equation (7)» This procedure 
"Will be eng)lo3nsd for the following "two cases; (l) slender cones for 
which 5, "the angle of flow inclination, is small congiared to 1 radian 
throvi^out the flow field, sind (2) cones for "which ci)-5, the difference 
be"tween the ray angle and the Inclination angle, is small compared to 
1 radian throughout the conical part of the flow field.^ 

Analysis of cones for which S<?C1 (slender cones). - With the 
restriction imposed "ipon 6 in this s"tudy, equation (o) may be reduced 
to the form 


^In general, (o-8 "will be small for blunt cones; however, at Mach num- 
bers "very large compared to 1, u-6 "will be small for relatively 
slender cones as well. 
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5 cotu+ ^ [1-M^sln^u) (1-6 cot w)^] =0 
dw 


(9) 


■v/here M is considered constant. This relation is still nonlinear in 
both 6 and u , and is not readily amenable to exact algebraic solution. 
An approximate solution may be obtained, ho-wever, in the following man- 
ner. Near the surface of the cone, equation (9) reduces to the linear 
equation 3 


•Which has the solution 


6 


cot (i) + 



6 = 


sinw 


( 10 ) 


No'ff 6 cot u decreases rapidly ■with increasing u away from the cone 
S'urface; hence it is suggested that -without appreciable loss of accuracy, 
the solution for 6 (near the s'urface) given by this equation can be 
substituted into the coefficient of dS/du in equation (9). Performing 
this operation, equation (9) becomes 

6 cot 0 ) + ^ [1 - M^in^w (l-k^ esc u cot co)^] = 0 (U) 


which is 
acc-uracy 


linear in 6 and can be integrated to yield (to the order of 
of this analysis) 


Sc 


1 + 


1 + 


2M^(sin^(j- sin^) 

1 + (/i+ij-M^sin^Uj,"” 
2M^(sin^c»)^,- sin^u) 
1 - i/i+HM^sin^WQ 


1 


( 12 ) 


which satisfies the boundary condition at the s'urface of -the cone 
(i.e;, 6 = 6^ at u = w^). 

Since M was ass-umed cons'fcant throughout the conical part of the 
flow field, it can conveniently be taken as the Mach number just down- 
s-feream of -the shock. In this case, flow conditions at the shock 
(i.e., Ms, Ug, and 6g1 are obtained by the simultaneous solution of 
equation (12) and the oblique shock -wave equations 


3 It is clear that cot cj co-uld be replaced ■with l/co in this equation^ 
however, the trigonometric operator is re-tained for consis'fcency -with 
the rest of the analysis in which such operators must, in general, be 
re"tained (to is restricted to be small only near the surface of the 
cone) . 
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2 _ (7+l)^Mo'^sin^a3g- It-(Mo^siii^-l) (7Mo^sin2us+ l) 


Ms‘^= 


[27 Mo^sin^ojg- ( 7 -I) ] [( 7 - 1 ) Mo^sin^Wg + 2 ] 


5r = tan ' 


cot g)g (Mo^sln^dig -l) 
2:^ Mo^-(Mo^sln^a)g -l) 


Having determined the flow conditions at the shock, the variation of 6 
with Cl) throughout the entire flow field about a cone is known from 
eq^uation (12) . 

It now remains to more accurately determine the magnitude of the 
velocity throughout the conical part of the flow field. For this pur- 
pose, equation (7) must he employed in combination with an expression 
relating m and 6. Equation ( 12 ) provides such a relationship; how- 
ever, it is unnecessarily conqilicated for determining the small changes 
in V. A much simpler expression of acceptable accuracy is obtained by 
neglecting 6 cot cii in the coefficient of dS/dca appearing in equa- 
tion ( 9 ) thus 3rielding 

S3. - -S coth) - . 

dco i-M^in%) ^ 


6 = k^i^sc^— (16) 

In this equation, kg eind are determined by the reqiiirements that 

S — be at ( 0 — Ci)p 


6 = 6g at (i)=ci)s 

where the latter quantities are known from the previous analysis. Sub- 
stitution of equation (16) into equation (7) provides the following 
relation for Y; namely. 


■^As will be shown later, in the cases where this siniplifying assumption 
introduces significant error in the variation of V with u (viz., 
when sin^ becomes of the order of 1) the restilts of the subse- 
quent analysis of cones for which 0)-5 « 1 may be applied with good 
accurac 
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V 


-kadS 




-k/(l-M^) 


1 + , 




5^-k 2 (i-m 2) j 


6d6 


(IT) 


which can he Integrated to yield 


Zn V = -kgZn 


5+/i^-k2^(l-M^) 


+ ^ ^ Zn [62_k22(1^2)] + kg (18) 


where kg is determined hy the requirement that V = Vg 'vrtien 6 = 5g. 
Replacing the constants in this expression with their values given hy 
the imposed boundary conditions, there is then obtained the relation 




■r 

[ (8^-5g^)cot^ (0c+ (5c^ -62)cot2o:^ 
[_ (^c^ ■ ^s^) cot%)g 


8c 2 - 8g2 

2 ' ' 
cot lOg-COt Wg 


_ /£cifsl_ 

/cot^ojg-cot^ojg + /i8Q^-8g^)cot^ti)B ~2V oot^cac-cot^cog 


8/^Ot%-COt®Ws + /(52-8s2)cot2(0c +(8g2 -62)cot2c0g 


(19) 


Knowing the velocity, 
the relation 


M 


the Mach number may 



} 


1 



of course,’ be determined from 



(20) 


and the pressure coefficient anywhere in the flow field may be obtained 
with the aid of the expression 



( 21 ) 
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where 


and 


^ _ 27 Mq^ sin^Ws -(7-1) 

■ ■ 

O 7+1 


( 22 ) 




(23) 


The Mach numher and pressure distributions (as veil as the orientation 
of the conical shock) throioghout the flow field about a slender cone are 
nov known. Kie flow field about a cone for which w-8 is small conrpeLred 
to 1 will be considered next. 

Analysis of cones for which a)-5«l. - In this case, equation (8) 
may be reduced to the form 


1 - (u-8) cot 0) + ^ [1 - M^(o)-5)^] = 0 (24) 

where M is again considered constant. This nonlinear relation will be 
solved in a manner analago'us to that emcployed in the solution of equa- 
tion (9 ). Ulus, near the svirface of the cone, equation (24) reduces toS 



or 


5 + (1 ) = 26q 

Conibining equations (24) and (25) yields 


(25) 


dS 


2 (o) - (i)(>) coto) -1 
1 - 4M2 (u-Wg)^ 


(26) 


which can be integrated to yield (substituting in the boundary condition 
at the cone surface) 

^t should be noted that -1 is the exact value of d5/dco at the surface, 
so long as the shock is attached. Consequently, in this analysis and 
the slender cone analysis both the magnitude and rate of change of 5 
with (1) are satisfied at the surface. 
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5 


fr -1 

a tan to^j sec ^tOg 

j tan tiy-( to-Ug) 

4}.i^an^t0c-i r ~ 

l+2M(o)-(Oc) 

tanuc 


tan tOp+ M(i+^ tan (Og) 
4M^( 1 - 2 M tan Uc) 


tan Wc+M(^ 

[l-2M(cj-o)c)] 4M^(if2M tan Wc) 


(27) 


Again M is ctiosen as tlie Mach number just downstream of the shock, 
thus enabling M (or Ms), tOs, and 5s to he determined by solving equa- 
tions (13) > (1^), aud (27) simviltaneously. Knowing these q[uantities, 
the relation between 8 and w throughout the remainder of the flow 
field may he obtained with equation (27) . 

In order to determine the small variations in V throughout the 
conical flow field, it is sufficiently acctirate (since u - 5«1) to 
assume a linear variation of 5 with to, namely, 

5 = k^u + kg (28) 

where k^ and kg are fixed by the requirements that 

5 = 8(j- at to =: (Oq 


and ' 


5 = 6s at (0 = ojs 

Combining equations (28) and (7) there is then obtained the .relation® 


V 


tan 


(8c“5) 


f ^s*^s 

\8c“8s 



(29) 


®The retention of the tangent operator in this relation is, to the 
accuracy of this analysis, optional. 
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■whicli is easily integrated to jdLeld 



cos 


cos (t»is-8s) 

(6c-S) (Mg-6s )~ 

Sc-«s _ 


5c -5s 


( 30 ) 


from ■wb.ich the Mach number and pressure coefficient anywhere in the flow 
field can he obtained hy employing equations (20) through ( 23 ) • 

There is a well-defined restriction on the range of applicability 
of the results of this analysis. In particular, when Ms(Ug-oJc) > l/ 2 , 
an imaginary value of 83 is obtained from equation ( 27 ) • This result 
permits the establishment of a boundary in the (Mq,8c) plane, given by 
Mg(a)s-ti3c)=l/2, separating the ar*ea in which equation (27) applies from 
the area in which it obviously does not apply. This boundary is shown 
in figure 2,"^ and it is clear that there is a minimum value of 8c for 
6iny Mo, below which the cone solutions Just presented do not apply. In 
the area below this boundary, the slender cone solutions must be employed. 


Flow Downstream of the Vertex of a Body 


Simplified expressions have been obtained for calculating the flows 
about cones operating at high supersonic airspeeds, and as was pointed 
out previotisly, these expressions can be 
employed to determine the fluid properties 
at the vertices of pointed bodies of revo- 
lution other than cones. The investiga- 
tion of the flow downstream of the vertices 
of such bodies is now undertaken. For 
this purpose it is convenient to first 
obtain an exact expression describing the 
fluid motion along a typical streamline in 
axially symmetric flow. To this end, con- 
sider sketch (a) showing the first and 
second family characteristic lines passing 
through a point on such a streamline . The compatibility equations defin- 
ing the variation of fluid properties along these characteristic or Mach 



'^The boundary curve was obtained by using the results of reference 8. If 
the equations Just developed were used, a boundary curve slightly below 
the one shown would have been obtained. Consequently, the indicated 
area of applicability of eqiiation (27) is slightly conservative. 
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lines may tie -written (see, e.g., reference l) 


-cos^p ^5 sin p sin 5 1 /dE'. 

- — = - — + + — ( — ) sin®p 

2^ + sin^p oci dci y • yR 


(31) 


along and 


cosS p ^ ^ ^ ^ 

Sco Scg y 7 RVdn 


+ sin^p ^"2 


sin^P 


( 32 ) 


along Cg. Comb inin g these eq;uations, and noting that 


and 



3p 

3ci ^ 3^ 

3s 

3cj^ 

3s 

3c2 

35 _ 

36 

^ci , 


3s " 

3ci 

3s 3cg 


Ss 




where 


“Ss Ss 


2 cos p 


there is obtained the relation 
cos3 p 


2^ + sin2 p 


d p = cos p<3S - 


Bs sin p sin 6 
+ L 

■Sci y 


dS 


(33) 


It is e-yddent that the, trooible some q,uantity in this eq,Tiation is 

the rate of change of flow inclination along a first family Mach line. 

An insight into the beha-vlor of this qiaantity in the flow field about a 
body of revolution operating at high supersonic airspeeds may be ob-tained 
however, -when the -val-ue of the similarity parameter K (the ratio of the 
free -stream Mach number to the slenderness ratio) for the body is large 
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compared to 1 .® To illustrate, consider first tbe flow "between the 
shock wave the surface in the region of the vertex of the body . 

Since the flow is conical in the limit as the vertex is approached, it 
seems reasonable to assume that the flow remains predominantly conical 
in type (see, e.g., the streamline pattern shown in sketch (b) for some 
distance downstream of the vertex. 



conical -type flow 
in region of vertex 


With this assunption® we have approximately, in the region tinder con- 
sideration (see anal3rsis of cone for which a>- 6 «l). 


dflJ 


-1 


where 


sincodn 
^(0 « * ■' « > ■ '> 

y 

and n is measured normal to 
the ray making the angle to 
with the X axis as shown in 
sketch (c). However, dn may 
be related to dc^ by the 
following expression 



dn = dC]^ sin (p + &-to) 


®t/hen K is less than 1 the second-order theory may ordinarily be 
ea^loyed to analyze the flow field about a body. 

®It is tacitly assumed throughout this and subsequent anal yses that 7 
of the fluid is greater than I. At the extreme hypersonic speeds 
where, for example, 7 of the air flow downstream of the bow shock 
approaches 1 , the streamline pattern must, as pointed out by Busemann 
in reference ll-, have the shape of the body in the infinitesimally thin 
region between the shock and the surface. In this limiting case, 
Busemann* s analysis will apply. 
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hence ^ combining these equations yields (in the notation of eqmtion (33)) 

( 31 ^) 


BS _ -sin O + S ~ (ii) sin o) 
Sci y 


To the accuracy of this analysis^ however (see equation (25) )^ 


Cl) = 25f» - 5 


Introducing this relation into the above expression for and 

recalling that 2(6 q-6) = (i)- 6«L, there is then obtained 


55 _ -(sin p sin 5 - 2(Sc-S) sin (6-p) - lj.(6c-S)^ cos p cos S] 

Bcj^ y 


Consistent with the assumption that K is large compared to 1, this 
equation may be further simplified, since the last two terms on the 
right in the numerator are negligibly small conpared to the first term. 
In this case equation (35) reduces to the form 


5s _ - sin 3 sin 5 
^Ci y 


(36) 


With the aid of this expression, equation (33) readily reduces to the 
more tractable form 


cos^ p 
+ sin2 p 


dp = (36 


(37) 


relating the change in Mach number with flow inclination along a stream- 
line in the region of the vertex. Equation (37) is recognized, of 
course, as the differential equation for Prandtl-^eyer flow. 
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Consider now the flow downstream of the region of the vertex (just 
considered) where the slopes of the streamlines are relatively s mall ^^ and 
the ordinates are relatively large. Flow in this region is, in general, 
certainly not predominantly of the conical typh, and hence the previous 
analysis cannot he expected to apply; however, for the values of K 
under consideration, certain ohse2rvatlons can he made regarding the rela- 
tive order of magnitude of the tenns in equation (33)- For example, it 
is noted that as K is increased for a given hody (which is tantamount 

to increasing LIq) , the term sin p sin S decreases while the term 

y* 

cos Bd6 increases. As K takes on large values congoared to 1, the 
former term must become, in fact, small compared to the latter term. It 
is thus indicated that for flow of the type under consideration equa- • 
tion (33) may he siiiq)lified to the two-dimensional form 


cos3 3 08 

dp = cos pdS dS 

21^ + sin^ B 
2 

• Similarly, the compatahility equations reduce to their two -^dimensional 
forms, which implies, of course, that the flow field itself is approxi- 
mately two-dimensional. This being the case, howevei^ it may he sho'vm 
(from the results presented in reference 9) that Bb/oci is genei*ally small 
compared to and hence that (58/Sci)dS is likewise small com- 

pared to cos pdS. Thxis essentially two-dimensional (Prandtl -Meyer) 
flow is found to prevail downstream of the region of the vertex. This 
result, althou^ remarkahly simple, is not entirely surprising since it 
wovild he anticipated that three-dimensional effects would he reduced as 
the local Mach numbers become large and hence the local Mach angles 
become small. In any case, we have the result that when K is large 
compared to 1, the Mach nu^er and flow-inclination angle along a 
streamline in the flow field about a hody of revolution are related by 
the familiar isentropic expansion relation (obtained by integrating 
equation (37)). 



^°In this disc\ission it is assumed that d5/Ss<0 on the body surface. 
In cases where Bs/SsX) on the stcrface, it is not evident that the 
argument will apply. It is suggested, however, that in these cases 
the previous arguments conceraing flow in the region of the vertex 
will apply over a larger region downstream of the vertex. 
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where A and B are different points on the same streamline. 

If the streamline flow pattern were known throughout the flow field, 
then the Mach n\miber distribution coxild he readily obtained with equa- 
tion ( 38) . Perhaps the most Tisef ul application of this and previous 
results of the analysis is, however, to the determination of Mach nuniber 
and pressure distributions at the surface of a body. Attention is there- 
fore tu3med to this calculation. 


Flow at the Surface of a Body 

The procedure for determining flow conditions at the surface is 
entirely anal ogous to that employed in the application of the so-called 
shock-expansion method to airfoils, and hence might be called a conical- 
shock-expansion method. For example, the Mach number on the surface at 
the vertex is obtained with equations (13), (l^), and (20) in combi- 
nation with eqtiations (12) emd (19), or (27) and (30) depending on 
whether 5c or the values of Mq and Bg are such that u-B «1, 

respectively.^^ The variation of Mach number downstream of the vertex 
is then obtained with equation (38) . Knowing the Mach number distri- 
bution, the pressure distribution (in coefficient form) on the surface 
is readily obteilned with equatioias (2l) throu^ (23). 


Simplified Expressions for Slender Bodies 


In the case of slender bodies the above described calcinations 
become so sinple as to warrant special attention. This additional sim- 
plification arises from the fact that now, not only is u-B small 
throughout the flow field, but also both uj and 8 are small through- 
out the field. Thus, too, M is everywhere large compared to 1, and it 
results, as shown in appendix A, that explicit solutions can be obtained 
for the Mach number and pressure at any point on the surface of a body. 
These solutions may be summarized (from appendix A) as follows. The 
local Mach number at any point is given by the relation 

” ' 7~T\ 

1 - (% 8„)(l- 

^^It is clear, of course, that the results of reference 8 may also be 
■used for this purpose. It is preferable, in fact, to "use these 
res\ilts when interpolation is not necessary. 
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Mo‘ 


_1 + ^ (MqSn)^ 


1 + 7 (MqSn)' 




(Mo5jj)‘ 


and 


(Mq8jj)‘ 


/8a\^ _ 

(.%y> \ + Z±1 


(Mo8u)‘ 


> (iw) 


-J 


The presstnre coefficient at any point on the sttrface may he obtained 
from the esjpression 




m 


It is interesting to note that these expressions predict the ratios of 
local to free-stream Mach numbers, and local to free-stream static pres- 
sures to he the same at corresponding points on related bodies, provided 
the flov fields about these bodies are related by the same value of the 
hypersonic similarity parameter. These predictions are identical to 
those of the hypersonic similarity lav (reference 7) and conseq.uently 
they provide a necessary check on the validity of the assumptions under- 
lying the development of the simplified methods of this paper. It 
remains nov to evaluate the accuracy of these methods by conparing their 
predictions vith those of more exact theories. 


Conparison of Approximate and Exact Calculations of Flow 
About a Body of Revolution and Discussion of Eesvilts 


The flow fields about a nvanber of nonlifting cones operating at 
stipersonic airspeeds have been computed nimoericalJ.y by Kopal (refer- 
ence 8) following a procedure similar to that first enployed by Taylor 
and MaccoU. The results of these calculations provide an accurate check 
on the approximate solutions developed in this paper for cones . 
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In figin^es 3 fhrougTi 6 6ire shown the variations of 6 vith w for 
cones having semiapex angles of 5*^, 10°, 20° and it0°. The predictions 
of the slender -cone solutions are observed to he in good agreement with 
the results of reference 8 for hoth the 5° and 10° cones (figs. 3 anfl k) 
even at Mach numbers as low as 1.5, and in from fair to good agreement 
for the 20° cone. The agreement iiproves with increasing Mach number, 
as would he expected. At the higher Mach numbers, where the cone solu- 
tion for which (»)-6«l applies, S is accurately predicted as a func- 
tion of 0 ) for all four cones, the range of applicability (in terms of 
Mach number) being the largest, of course, for the ^° cone. It is 
Interesting to note that fair agreement is obtained for the latter cone 
down to a free -stream Mach number as low as 2, which is only slightly 
above the Mach number for shock detachment. 

The variations of velocity with ray angle throughout the flow 'fields 
about the above described cones are shown in figures 7 throu^ 10. The 
agreement between the cone solutions of this paper and the results of 
reference 8 for these variations parallels that found for the flow incli- 
nation angle as a function of the ray angle. In general the velocity 
variations are small, decreasing percentagewise with increasing Mach 
nimiber, as was asstuned in the analysis. It is encouraging to note, too, 
that either or both of the two ajjproximate solutions generally provide 
an accurate prediction of the variation of 5 and V with to throu^- 
out the flow field (e.g., figs. and 8). It appears that a rule of 
thimib for choosing the preferable solution is to enploy the cone solu- 
tion for which to-8«l whenever it is applicable (fig. 2). 

The press\are coefficients , on the surfaces of the four cones have 
been calculated tising the two approximate solutions, and the results of 
these calculations, along with the predictions of second -o 3 rder theory 
(taken from reference 3) aud the results obtained from refesrence 8 are 
presented in figiire 11. It is observed that the two approximate solu- 
tions overlap to predict pressure coefficients for the 5°, 10°, and 20° 
cones that are in good agreement with those of reference 8 at Mach num- 
bers from approximately 1.5 to Infinity. Conparable accuracy is obtained 
from the cone solutions for which u)-5« 1 for the 40° cone at Ifech 
numbers above about, 3. In general, of course, the agreement in^iroves 
with increasing Mach number. The second -order theory yields more accu- 
rate results than the slender cone theory for the 5° 10° cones at 

the, lower Mach numbers, althou^ the reverse seems to be the case for 
the 20° cone. Neither second -order nor slender -cone theory is appli- 
cable to the ^° cone. 

From the preceding conparison of the conical flow' calculations of 
reference 8 with the predictions of the simplified solutions of this 
paper, it is indicated that the latter solutions may be enployed to 
predict the properties in the flow field about a cone with from good to 
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excellent acctnracy, depending on -whether the Mach number of the free 
s-fcream corresponds to intermediate or high supersonic airspeeds. It is 
therefore suggested that these solutions may, for example, he particu- . 
larly useful at high siiper sonic speeds for accurately, determining the 
conical flow fields about cones at Kfeich numbers not trea-ted in the M.I.T. 
■tables (reference 8) . 

It remains now to determine the accuracy with -which the solutions 
for flow about cones in combination -with the isentropic expansion equa- 
tions predict the flow at the surface of bodies of revolution other than 
cones. The Mach number and pressure distributions at the surface of a 
family of ogives operating at a free -stream Mach number of 6 were calcu- 
la-ted using the simplified methods of this paper in the manner described 
in the analysis. These distributions are presented in figures 12 and 
13 for -values of K vaiying from l/2 to 2. Also shown are the results 
presented ir. reference 10 , ob-bained. -with the method of characteristics 
including effects of ro-fcation. In general it is observed that -with 
increasing K, the conical -shock -expansion theory of this paper yields 
Mach numbers and pressures that are in better agreement -with the predic- 
tions of charac-fceristics theory. Indeed at K as low as 2, the pre- 
dictions of the two theories agree -within the accuracy of the character- 
istics solution on the ogi-val nose.13 The several asstimptions made in 
■the development of the conical -shock -expansion method are thus Justified 
in part. It is somewhat surprising, however, that the method works as 
■well as it does at -valties of K near 1. It is suggested that this 
res\xlt may be pec-uliar to streamlines near the surface of a body, for it 
may easily be shown (with an argument analogous to that presented in the 
analysis) that in the region of the ver-tex -where three-dimensional 
effects are perhaps most pronounced, flow in stream -tubes adjacent to 
the surface nevertheless shows agreement -with the -two -dimensional 
(Prandtl -Meyer) relation be-tween stream inclination and Mach number or 
pressxire for -val-ues of K approaching 1. 

The slender body solution also displays increasing accuracy -with 
increasing K, although it consis-tently predicts -too high Mach numbers 
and too high press-ures on the sirface near the vertex. Ihis result is 


^^The Mach number dis-tributions presented here were calculated from the 
pressure dis-tributions given in reference 10. 

^3The frequently sugges-ted method of determining surface pressures on a 
body of revolution by as sinning the pressure at a point -vrLU be the 
same as on a cone -tangent to the body at that point -was eilfeo -tried 
and, in general, the press-ures were too high. The error was less at 
larger K, howe-ver, although it -was greater than for the conical- 
shock -expansion theory. As -would be expec-ted, the Mach number distri- 
butions -were generally in considerable error. 
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primarily a consequence of the approximate nature of the conical flow 
solution en 5 )loyed in the development of the theory. In any case even 
at a K of 1, the slender hody theory displays sufficient accuracy for 
many engineering purposes. Indeed, upon inspection of the pressure dis- 
tributions, it is evident that, ailthough doubtless fortuitously, the 
slender body theory will yield more accurate drag coefficients than the 
conical -shock -expansion theory at the lower values of K. 

Comparisons similar to those just discussed were made for the other 
ogives considered in reference 10, and in general the same results were 
obtained; namely, for K>1, the simplified thebries of this paper were 
in good agreement with the predictions of the characteristics solutions, 
the difference between the predictions of the latter theory and the 
conical -shock -expansion theory being of the order of accuracy of the 
characteristics solutions at K as low as 2. 


.CONCLUDING EEMARECS 


With the assumption that the flow at the vertex of a body is coni- 
cal, it was found that simple approximate solutions can be obtained 
which yield the Mach number and pressure at the surface over a consider- 
able range of free -stream Mach numbers and apex angles. In the special 
case of cones, these solutions define the entire flow field with good 
accuracy, and may therefore provide a useful adj-undt to the well-known 
M.I.T. tables for flow about cones. 

The investigation of flow downstream of the vertex revealed that 
when the hypersonic similarity parameter K (the ratio of free -stream 
Mach number to slenderness ratio) is large compared to 1, the Mach num- 
ber varies with flow inclination angle along a streamline in approxi- 
mately the same manner as for two-dimensional (Prandtl -Meyer) flow. In 
the special case of streamlines near the surface, it was suggested that 
this parameter may approach 1. This and preceding results obtained for 
flow at the vertex were combined to yield what might be termed a conical 
shock -expansion method for calculating the Mach number and pressure 
distributions at the surface of a body. In the special case of slender 
bodies, exceedingly simple explicit expressions were obtained for these 
quantities. 

Surface Mach number and pressure distributions were calculated with 
the simplified methods of this paper for a family of ogives traveling at 
a Mach number of 6. These distributions were in good agreement with 
those obtained with the method of characteristics at values of K 
greater than 1. In the case of the conical -shock-expansion calcula- 
tions, the agreement was within the order of accuracy of the character- 
istics solutions at K as low as 2. 
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Because of the relative simplicity of the proposed methods of deter 
mining the flow at the surface of a body operating at high supersonic 
airspeeds (K >l) , these methods should prove useful for engineering pur- 
poses. It is also suggested that the same general approach may he 
applicable (again for K >l) in developing methods for calculating flow 
at the surface of pointed bodies of revolution at small angles of attack 
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APPE3TOIX A 

FLOW AT THE SUEFACE OF A SLENDER BODY OPERATING 
AT HIGH SUPERSONIC AIRSPEEDS 

If a slender "body (i.e., a body on tHe surface of ■wLic'h the slopes 
are everywhere small con^jared to l) is operating at free -stream Mach 
numbers very large compared to 1 (again, of co^rrse, K>l), the local 
Mach numbers will likewise be large compared to 1. It follows, then, 
that the inclination of the nose shock wave will be small and, conse- 
quently, that w will always be small. In this case the relation 
between 5 and u at the vertex is extremely simple; namely, (see 
equation (lO)) 



Combining this expression with equation (7), the relation 

Vs 'VS/ 


(Al) 


(A2) 


defining the velocities in the flow field is easily obtained. Hence, 
the Mach number on the s\irface at the vertex, %, may (to the order of 
accuracy of this analysis) be related to Mg by combi n i n g equation (20) 
with this expression to yield 



Now the oblique shock -wave equations for flow of the type under 
consideration reduce to 


The equations developed in the following analysis could be obtained 
from the final expressions previously developed by reducing them 
to conform to the assumptions of this slender -body theory. It seems 
sinpler, however, to proceed, as indicated, from the basic -flow 
equations . 
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(7+l)%^(Motos)' 


M ^ = 
^ s 


[27 (Mows) - (7-1)3 [(7-l)(MoC0s) + 2] 




and 


Pq 7+1 7+1 


(A?) 


vtiere 


Mo.s = Ws (2^ |l . JziMlii I 

^ t (7 +i)Mo5e 


(a6) 


Combining equations (Al) and (a 6) there then results 


MotOs 




(MoS]j)‘ 


(AT) 


and 


Ss 

5n. 


(Mo5n)' 


1 + •2~(Mo5jj)^ 


(a8) 


Equation (A^) can now he written (considering equation (A7)) 



Mq2 

1 + 



"l + 7 (Mo5n)^ 


1 + (^^w)^ 


(A9) 


Equations (A3), (a 8), and (A9) provide the Mach number on the sur- 
face at the vertex. IQiowing 1%, a simplified expression for the Mach 
number anywhere on the svirface of the body may be obtained from equa- 
tion (37)* Since the loceil Ifech nunlier is assumed very large, this 
equation may be reduced to the form 


2 


dp = d5 
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wTiicli, in turn, can "be written, upon integration, to yield the local 
Mach numher 


M = 



Now the pressure coefficient is given hy the expression 



1 . 

\5o Ps 



(AlO) 


(All) 


The pressure rise across the shock may he obtained hy combining eq.\ia- 
tions (A5) and (A?) to yield 

^ = 1 + 7 (Mo5uf 

Similarly, the ratio of the pressure at the surface to the pressure at 
the shock can he expressed (to the order of accviracy of this analysis) 
in the form 



VJith the aid of these expressions, equation (All) can now he written 


P = 


2 

7Mo^ 


1 + 7 (MqBn)^ 


-BZ 

7-1 



yielding the pressure coefficient at any point on the surface of the 
body. 


(A12) 
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Free^stream Mach number. Mo 


Figure 2.-Variatlon of cone half-angle with free-stream Mach 
number showing range of applicability of cone solution for 
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(a) Mo=L5058 



Ray angle,v,degrees 
(c) Mo=l3.l96 



Figure 3 .—Variation of flow inclination angle with ray angle for 
a 5* half-angle cone operating at Mach numbers of 1.5058, 
4.8602, and 13.196. 
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Figure 4 .—Variation of flow Inclination angle with ray angle for 
a !0° half-angle cone operating at Mach numbers of t.5957, 
5.4223. and 15.146. 
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Figure 6 .—Variation of flow inclination angie with ray angle for 
a 40“ half-angle cone operating at Mach numbers of 2.0242, 
5.8336. and 13.101. 




MCA TTT 2579 



(b) M<,=4.8602 


.9810 

.9808 

.9806 

.9804 


.9802 


1 1 

Equation (19) 

Eauation (30) 






/ 

?eferenc 

'e 8 



- 






























5.50 575 6.00 6£5 

Ray angle,o),c/egree$ 
(c) Mo>^l3.l96 


6.50 6.75 7.00 


Figure 7 .—Variation of the ratio of local velocity to limiting 
velocity with ray angle for a 5® half -angle cone operating at 
Mach numbers of 1.5058, 4.8602, and 13.196. 
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Ray angie,o>, degrees 
(c) Mo =15.146 

Figure 8 . — Variation of the ratio of iocai veiocity to timittng 
vetocity with ray angte for a tO* haif-angte cone operating at 
Mach numbers of t.5957, 5.4223. and t5.t46. 
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Figure 9 —Variation of the ratio of local velocity to limiting 
velocity with ray angle for a 20° half -angle cone operating 
at Mach numbers of 1.4672, 5.5457, and 9.5928. 
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(b)Mo^5j8336 



Ray angle, w,degrees 
(c) Mo =13.101 



Figure !0 .—Variation of the ratio of local velocity to limiting 
velocity with ray angle for a 40^ half -angle cone operating 
at Mach numbers of 2.0242, 5.8336, and 13.101. 
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Figure 12 .—Variation of local Mach number along ogives 
corresponding to four values of the similarity parameter, 
K, for a free— stream Mach number of 6. 
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(a) K^O.5 
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Longitudinal coordinate, percent nose length 
(b) K = i 



Figure !3 —Variation of pressure coefficient along ogives corresponding to four 
values of the similarity parameter, K, for a freest ream Mach number of 6. 
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Figure 13— Concluded. 
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